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^ ■ Abstract 
(N , 

. The thermal evolution of the spectral densities derivable from the two-point 

i functions of the elementary and the quadratic composite fields of the 0{N) model 

' is studied in the isosinglet channel and in the broken symmetry phase at infinite 

r"| ' N. The results are applied with realistic parameter values to the = 4 case. They 

Q-i! provide a reasonable description of the a meson at T = 0. Threshold enhancement 

Q^' is observed around T ~ l.OTmTr- For higher temperatures the maximum of the 

^ ■ spectral function in the single meson channel decreases and becomes increasingly 

• • , rounded. 

X 

^ ; 1 Introduction 

There is increasing interest in finding signatures for the characterisation of strongly in- 
teracting matter around the transition temperature from the hadronic to the quark-gluon 
regime. Density and temperature induced features of the correlations in the soft pion 
spectra is a frequently discussed issue. The majority of these pions comes from resonance 
decay. One expects that temperature induced shifts in the position and width of the p 
and a pole location is sensitively reflected in the spectra. 

The effects related to the variation of the location and the width of the broad a 
resonance were investigated in the past decade systematically The main theoretical 
framework of these investigations was the 0(4) linear sigma- model. Various reorganisa- 
tions of the perturbative series were proposed for the self-energy calculations. Despite 
of the rather large value of the nonlinear coupling A^j ~ (9(30 — 70), one loop calcula- 
tions gave important qualitative hints for the enhancement of the spectral function in 
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the (T-channel near the two-pion threshold as the temperature approaches the transition 
region. 

In the non-relativistic context it has been recognised long time ago that in scalar 
theories the hybridisation of the two-point functions of the elementary and quadratic 
composite fields play essential role in the determination of the elementary excitations in 
the whole temperature interval of the broken symmetry phase. Neutron scattering data on 
superfiuid helium were successfully interpreted by taking into account the hybridisation 

Large techniques are known to account for this phenomenon correctly, while ordi- 
nary perturbative approaches usually miss it 

An additional point in favour of the use of the large N expansion is that it leads to 
a second order transition point in agreement with the well-known behaviour of the 0{N) 
modeip]. The lowest order computations done with any (even optimised) version of the 
conventional perturbation theory yield a first order phase transition when the temperature 
is increased. 

Large techniques are extensively studied also in connection with the out-of- 
equilibrium evolution and thermalisation of relativistic Bose-condensates 0, ||, ^, 0, [n|] , 
though the hybridisation in this respect did not receive till now much attention. 

Our main goal with this investigation was to apply large A^ techniques for finding 
the finite temperature variation of the spectral function in the elementary and quadratic 
composite isosinglet channels, which influences the vr — vr correlations measured in heavy 
ion collisions to a large extent. For this purpose the parametrisation of the linear 0(4)- 
symmetric meson model is chosen to reproduce as closely as possible the accepted zero 
temperature fit to the complex mass, M„ — iV /2 of this broad resonance fT^ . 

In this Letter, after summarising the general formalism, mostly the results of the 
chirally invariant limiting case will be discussed, where the pions are massless, since the 
main ideas can be presented in this case very transparently. Results relevant for the case 
with explicit symmetry breaking will be presented shortly in the concluding part of the 
paper. 

We shall work in the leading order large A^ approximation to the Schwinger- Dyson 
equations of the 0(4) model in the broken symmetry phase. The effect of next-to- leading 
order corrections is the subject of our ongoing research. 

2 The model and its renormalised dynamical equa- 
tions 

The Lagrangian density including a term reflecting explicit breaking of the 0{N) sym- 
metry is the following: 

L = lid.rd^r - m'rr] - ^i<Pi\<p')' + ^h'. (i) 

The broken symmetry phase can be studied after an appropriate shift, by introducing the 
symmetry breaking background: 

(/.'^ ^ (ViV$ + 0\0O. (2) 

The quantum field representing the fluctuations of the order parameter is termed o", while 
the modes transversal to it are the Goldstone pions. 
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The equation of state, which determines the absolute value of the condensate comes 
from the requirement of the vanishing quantum expectation for the coefficient of the term 
linear in 0^ in the shifted Lagrangian: 
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A w,„^0^ h 



The quadratic fluctuations of the shifted flelds are computed with an accuracy 0{N^), 
therefore only the contribution of the pions is retained. Anticipating a non-zero Goldstone 
mass due to the explicit symmetry breaking (which will be verifled below) one flnds the 
following relation connecting renormalised quantities (only contributions proportional to 
are kept): 



m„ -\ <P{I ) -\ -m}.{l)m — ^^^-^ / ay— = — ; — -. 4 

^ 6 ^ ' 967r2 ^ M2 127r2 J^np. ev - I ^T) ^ ' 



Here the following renormalised couplings were introduced (momentum cut-off was applied 
for the regularisation of some divergent integrals): 

m% eA^ m% 1 1 , eA^ 1 

I £i_ln = a I In = 

X ^ 967r2 967T^ Mi Ar ' A ^ 967r2 A^' ^ ' 

The issue of the normalisation point Mq will be discussed below. 

At leading order [N = oo) the only contribution to the self-energy of the Goldstone 
modes comes from the tadpole diagram. Therefore one has: 

Ga\p) =p'-m^- ^$^(T) - ^{{rr) = - (6) 

Here the sum representing the mass term was simplifled in view of the equation of state 
(|). This equation implies that at A^ = cxd the Goldstone particle is stable and its mass, 
mg(T) = h/^{T) increases with the temperature. 

With this identiflcation the renormalised equation of state (^) can be cast into a more 
practical form by eliminating and h in favour of the T = value of the condensate 
($o) and of the pion mass {itigq): 



+ r, dy^f^^'iey - ir = ^ - 1. (7) 

This form makes it clear that after Xr and Mo/mco are chosen in the process of renor- 
malisation at T = 0, one flnds from this equation mciT) /itlgo as a function of T/itlgq if 
the physical input ^^/h = f^/Nml is made. 

The cr propagator receives C(A^°) contribut ion from the inflnite iteration of the bubble 
diagrams b{p), where on both lines Goldstone flelds are propagating: 
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Without entering its derivation we give here also the expression for the leading large N 
propagator of the quadratic composite field — ((0")^))(x, t): 

^ {p'-hmT))b{p)/6 + <l>\T)/3 

(p2 - /./$(T))(1 - \Rb{p)/Q) - A^$2(T)/3- ^ ^ 

It has the same denominator as Gh{p) making explicit the hybridisation of the two objects. 

The bubble contribution with external momentum po, p is the sum of a zero temper- 
ature and a T-dependent part, b{p) = h^ij)) + hripo, p), which read as follows: 
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-{rii - n2) 



Pq - uji - UJ2 + It pq + uoi + UJ2 + te 
1 1 
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_Po - iOi + UJ2 + te po + uji - 002 + «e. 

where = 1 / {exp{(3uJi) — 1) and Ui = (q^ + m^)^/^,cj2 = ((q + p)^ + m^)^/^. The first 
term in the expression of bo{p) is cancelled in the expression of the a propagator by the 
divergence of the bare coupling A and the inverse propagator, expressed in terms of the 
renormalised quantities is finite. Note, that bo{p) has to be evaluated with mG{T)\ 

The phenomenologically most interesting object is the spectral function of the order 
parameter field a, defined as 

Ph{po,P,T) = --lmGHiPo,P,T). (12) 

71 

For non-zero T pn has a singularity at po = due to the Bose-Einstein factor n{(3po/2). 
In the chiral limit h = this point coincides with the two-pion threshold. In order to 
make the effects related to the physical excitations visible around this point in Fig.|I| we 
show pi(po) 0, T)T^ = (1 — exp(— po/2T))p/^(po, 0, T)T^. The curves are drawn for various 
temperatures below and a representative value of Xr = 310. 

The temperature dependence of the order parameter needed for its evaluation is ob- 
tained from the equation of state which simplifies in this case Q to 

*^<^'-ifl-^). T^.UK (13) 



T2 12 V T2 

The shape of the spectral density starts with a broad bump at T = which is shifted 
towards lower frequencies when the temperature is increased. Already for T = pn has 
finite value at po = 0, where its value is gradually increasing. At a temperature T ~ O.TT^ 
the bumpy structure at finite po completely disappears and a simple Lorentzian shape 
develops with its maximum located at the threshold. 

This evolution can be transparently interpreted if the temperature dependent location 
of the physical pole of the a propagator is found. The physical pole is located in the lower 
half of the complex pQ-p\a.ne, therefore one is faced with the analytical continuation of 
bxip) into the lower po-halfplane, since originally it is defined only in the upper halfplane 
by the Landau-prescription: e > 0. As becomes clear, this threshold enhancement is the 
manifestation of the gradual chiral symmetry restoration as T tends to Tc. 
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Figure 1: Temperature dependence of the modified spectral function pi in the a meson 
channel for h = 

3 The physical quasiparticle excitation in the a chan- 
nel 

Below we give some details concerning the temperature variation of the excitation spectra 
and the interpretation of the spectral function with its help in the chiral limit {h = 0). 
One starts with the analysis at T = 0, where one fixes the renormalised parameters 
m^, Aft, remaining after h was set to zero. The renormalised mass is related to the size of 
the condensate (or in view of ([T3| ) to Tc) by the well-known equation: — m|./$o = ^r/^- 
Xji will be chosen from a range where one finds for T /M„ values which are the closest 
to the experimental one. For this we determine the T = pole of the a propagator by 
looking for the zeros of in terms of renormalised quantities: 



A 



P 



R 



1 



ARln(-j9VM2)/967r2 



0. 



(14) 



The physical solution at rest of this equation is parametrised by putting p = 0,po = 
Mq exp(— ii^o)) 7'"/2 > v^o > 0. (One can find its mirror solution in the third quarter). The 
notation shows, that the absolute value of the pole is chosen for the normalisation point 
defined in (^ and used in (0). If a physically satisfactory solution is found for some value 



Aij with this normalisation point, the renormalisation group invariance of (|Tj) ensures that 
for a different choice of Mq the same ratio T /Ma- is found at some appropriately shifted 
value of \r. In this way the actual value of \r cannot be said to be large or small! 

Since in view of (^) in the chiral limit one has a very simple equation for the critical 
temperature, our solution provides the mass = Mq cos ipo and the imaginary part 
Mq sin (fio in proportion of the critical temperature. It is interesting to note that using for 
2$ the experimental value of f-„, eq.(|13D gives Tc = 161 MeV, while the lattice 
'3 ± 8) MeV [1T3. The agreement is much better than expected. The 
real and imaginary parts of the complex physical pole are shown in Fig.|^ as a function of 
Xr. 



simulations yield 
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Figure 2: The imaginary and real parts of the physical poles at T = in the chiral limit 
and also for h ^ 0. Also shown is the logarithm of the tachyon pole position in proportion 
to the mass of a for various temperatures below Tc and for h = 



In addition to the physical zeros of the inverse propagator scalar theories are known 
to have a tachyonic zero on the positive imaginary axis (po = iM^) (see [§ and references 
therein), which makes the theory unstable. Therefore the theory can be used as an 
effective theory until the physical scale Mq is acceptably smaller than Ml. The position 
of the tachyonic zero is also found from ([T^) . (The fact that the tachyonic pole shows 
up also in Gh not only in F is a manifestation of hybridisation). It is clear from Fig.^ 
that for A/j < 400 one has Ml/Mq > 4 and the effective approach is well justified. The 
representative value \r = 310 corresponds to a ratio M„/T = 1 at T = 0. The value of 
is estimated to be 7<l>o ~ 350 MeV. 

At finite temperature we are interested in the temperature dependence of the cx-pole 
at rest, therefore we set p = in (|11]). The analytic continuation onto the lower plane 
can be carried out by adding to the expression of fey valid in the upper halfplane a term 
which makes its imaginary part continuous when one approaches the real po-axis either 
from the upper or the lower halfplane: 



n — 1 



X Z + X 



— , Im po > 0, 



i 1 

briPo) = b^iPo) - 7, Im Po < 0, Re po > 0. (15) 

47r exp(z) — 1 

The roots are parametrised the same way as in the T = case, pq = M{T) exp{—i(f). 
After one determines $(T) from Eq.(0), from the complex equation 

1 

^ ^247rexp(Mexp(-i^)/2T) - 1 ~ ^ ^ 

one can find the complex pole as a function of the temperature. All quantities are mea- 
sured in units of ~ $0- 
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Also at finite temperature one can study the tachyon solution in the upper halfplane 
using b^liM^). One might wonder if the temperature dependence of the tachyonic pole 
does not restrict further the coupling range where the effective use of the scalar theory 
is consistent. The broadening in Fig.|] of the line corresponding to the tachyonic pole 
reflects the slight decrease of its mass scale with increasing temperature. However, it is 
clear that in the region A/j < 400, where the effective approach is consistent for T = 0, 
the position of the tachyon pole is practically unchanged. 

Similarly, one can look directly for physical roots along the negative imaginary axis. 
Then using from (^) one finds that above a well-defined Timagi^n) < the physical 
pole becomes purely imaginary: —iM{T). If fi = M{T)/2T << 1, for the highest 
occupied low frequency region, one can use the expansion of the Bose-Einstein factors 
appearing in the last two terms of ([T6|) with respect to the powers of fi and keep the 
leading terms: 

487r2 Mo 3 AT^ /i2 487r/i ' ^ ^ 

This equation could have been derived directly, if one would have applied from the be- 
ginning the same approximation in Eq.(|rT]). Clearly, the region around the origin (the 
critical region) can be analysed also directly with help of this simpler equation. Note 
that one has to go one step beyond the classical approximation to achieve a consistent 
approach. 

The general pattern of the trajectory of the physical pole with increasing temperature 
was the following. The real part of its position started to diminish when T was raised. 
In this way the broad bump of the spectral function moves towards the origin, and its 
width is increasing due to the slight increase of its imaginary part. Depending on A^j the 
root reaches at some Timag < Tc the negative imaginary axis and collides with its mirror 
root arriving from the third quarter. Each one is joining smoothly the trajectory of one 
of the pair of imaginary roots which appear just at Ti^ag and move opposite directions. 
We observe that the a-bump gets lost in the background before the temperature reaches 
Timag- The root approaching the origin along the imaginary axis produces in the spectral 
function a shrinking shape which becomes a Lorentzian only very close to Tc. Eventually 
for T = Tc it builds up a term proportional to 5(po)/po- 

We see that the complex pole of the cx-propagator qualitatively accounts for the be- 
haviour of the spectral function near and below the critical point as well. Note that the 
other pole moving away from the origin along the imaginary axis as well as the tachyon 
disappear from Gh at T^. They remain poles of F only. 

In the vicinity of the critical point, where the ^ = T/87r$^(T) is the dominant length 
scale and the condition Po? |p| << T is fulfilled, one can derive an equation also for the 
soft modes with nonzero momentum: 

-3.|pk?i^ln^-^M^,„?i^.l. ,18) 

|p|2 Po + IpI Atc Tc T^ 

Its solution in the approximation, when on the left hand side only the first term is retained 
exhibits the form of dynamical scaling [|14|, |T^: po = |p|^/(|p|0 with z = 1. The second 



term provides the leading correction to scaling. In 0{N) models, the dynamical exponent 
z = d/2 has been obtained for finite on the basis of scaling and renormalisation group 



arguments, where in our case d = 3 [jT6|, 0. For the correct interpretation of the 
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Figure 3: The spectral functions for the hnear sigma model with realistic couplings and 
explicit symmetry breaking. On the left figure the imaginary part in the a-channel, on 
the right the same for the composite channel is shown 



situation it is important that there are two distinct hydrodynamical regions in the 0{N) 
model for large . (The system studied in |jT6| can be regarded as a lattice regularised 

version of the linear a- model). In the true critical region z = d/2 is valid, in a precritical 
region z = 1 — 8Sd/Nd + 0(l/iV^), where 5*^ = 2/7r^ for d = 3. The first region shrinks 
when becomes large and completely disappears at iV = oo. Furthermore, it has been 
found that z agrees with the dynamical exponent determined for the quasiparticles, that 
is outside the hydrodynamical regime, at least to ^^(l/A^) [|l9l. Details of the application 
of this analysis to the present case will be published in our forthcoming paper. 



4 The case of the expHcit breaking of chiral symme- 
try 

Finally, we wish to discuss shortly the results of the application of the leading large 
solution of the 0{N) model to the low lying effective n — a system where an explicit 
symmetry breaking is necessary. 

Since in this case one cannot speak of any phase transition, it is more sensible to look 
for the roots of G]j^{po, 0) in proportion to $0 = /7r/2. The value m^Q/$Q ^ 9.06 is fixed 
by phenomenology and choosing A^j with the restriction Mq/Mi < 4, one can find the 
zero temperature position of the cr-pole. In Fig.^ its real and imaginary parts are already 
shown as a function of Xr together with the results obtained for h = 0. In the region 
of Xji allowed the ratio M^/T^ moves away from the phenomenologically preferred range 
(Mo- = 3.95/^, Mo-/r ~ 1.4) emphasising the need for a next-to-leading order calculation. 

Now $(T) is determined from (^. The spectral function of the single a-channel is 
shown in the left of Fig.|^ for various values of the temperature in the transition range 
T = (0 — 1.2)mGo with = 400. The location of its maximum at T = is correlated 
with the location of the T = pole. Maximal threshold enhancement is observed in 
the close neighbourhood of 1.074mGo with the spectral function numerically exhibiting 
a cuspy form, with a very narrow width. However, when the temperature is further 
increased, the maximum is pulled away from the actual position of the threshold and a 
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broad rounded structure is seen. We return to the detailed interpretation of this structure 
in terms of the pole trajectories in a forthcoming publication. 

As we have mentioned above also the spectral function belonging to the propagator 
(^) of the composite field can be constructed to leading order in N. The corresponding 
spectral function reflects the excitations in the composite (0")^ channel. Its temperature 
dependent deformation is displayed on the right hand of Fig.^ for the same (realistic) 
parameters as used in case of pn above. Its behaviour is rather similar to the single 
meson spectral function. No signal of any meson-meson bound state can be observed. 
The only difference one notices is that for large frequencies this function approaches a 
constant, which reflects mostly the large residuum of the tachyon pole in this propagator. 

5 Conclusions 

We have described the temperature dependence of the elementary excitations of the 0(4) 
model in the leading order of the expansion. For the chirally symmetric case a very 
suggestive picture of the complex pole evolution makes unique the interpretation of the 
change of shape of the single-particle spectral function when the temperature approaches 
the critical point. 

The restrictions on the range of the scalar self-coupling Xn arising from the requirement 
of keeping distance from the tachyonic pole, for a fixed normalisation point Mq prevent 
us in finding a fully realistic a-particle mass at T = when the physical pion mass is the 
input. We have found that the threshold enhancement in both the elementary and the 
composite spectral functions is maximal at some Tenh < ^go- Beyond this temperature 
the cuspy maximum becomes rounded again. We believe that this qualitative feature 
remains valid when the next to leading order corrections will be included. 
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